Accretion of phantom scalar field into a black hole 
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Using numerical methods we present the first full nonlinear study of phantom scalar field accreted 
into a black hole. We study different initial configurations and find that the accretion of the field into 
the black hole can reduce its area down to 50 percent within time scales of the order of few masses 
of the initial horizon. The analysis includes the cases where the total energy of the space-time is 
positive or negative. The confirmation of this effect in full nonlinear general relativity implies that 
the accretion of exotic matter could be considered an evaporation process. We speculate that if 
this sort of exotic matter has some cosmological significance, this black hole area reduction process 
might have played a crucial role in black hole formation and population. 

PACS numbers: 04.70.-s,95.36.+x,98.80.-k,04.25.D- 

The study of black holes in general relativity as ther- 
modynamical systems is an extensive area of research. 
One of the most interesting features of black holes that 
can be studied, is the behavior of the area of the hori- 
zon during the process of accretion of matter. Assuming 
that the null energy condition is satisfied, the area of the 
horizon can only increase [l[ . If this condition is violated 
then "strange" things can happen, for example, the area 
of the horizon could decrease and the black hole could 
shrink. This scenario could be astrophysically and cos- 
mologically interesting if some evidences of exotic matter, 
violating the null energy condition would be found. For 
instance, nowadays standard cosmological models, based 
on observations allow that matter with these types of 
properties might well play the role of dark energy Q . 

In this article we present the accretion of a phantom 
scalar field into a black hole in the full nonlinear case 
using numerical methods. The model of phantom scalar 
field is based on the following Lagrangian density 



implications in astrophysical scenarios, because the area 
increasing theorem does not apply in this case 

Previous experience shows that such a scalar field with 
V = is able to produce an area reduction effect, as 
found during the evolution of wormholes supported by a 
scalar field that collapse to form black holes [3( ■ Inspired 
by such an area reduction process we explore in this arti- 
cle the full nonlinear accretion of a phantom scalar field, 
that is, the case where a nonzero potential is involved. 

The evolution equations. We formulate Einstein's 
equations coupled to the phantom scalar field in a suit- 
able form for the numerical integration of spherically 
symmetric black hole solutions. We use geometrized 
units for which the speed of light and Newton's constant 
are equal to one. We write the metric in the general form 
for the coordinate system (t, r, 9, ip) in spherical coordi- 
nates: 



(i) 



where 1Z is the Ricci scalar of the space-time, c/ M „ is the 
space-time metric, <j) is the scalar field and V(cj)) its po- 
tential, which we choose to be V(4>) = Xcj) 2 with A a given 
constant. The difference between the phantom and the 
standard scalar field is the sign in the kinetic term in the 
Lagrangian. 

When the action constructed with such a Lagrangian 
density is varied with respect to the metric, the aris- 
ing Einstein's equations are related to the stress energy- 
tensor Tfj, v = -d^dycj) - ^g f _ lI y[-d a cl)d a (j> + 2V]. Such 
a stress-energy tensor violates the null energy condition 
T^ v k^k v < 0, where k^ is a null vector. The immedi- 
ate implication of this violation is the violation also of 
the weak energy condition, which in turn implies that 
observers following timelike trajectories might measure 
negative energy densities. 

The fact that a phantom scalar field violates the null 
energy condition motivates the study of possible unusual 



ds z = -[a'- (3 r (3 r »— dt 2 + 2(3 r ^dtdr 
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+ - \g rr dr 2 + g ee {d6 2 + sin 2 Odtp 2 )] , (2) 
X 

where x acts as a conformal factor relating this metric to 
a spacelike flat metric, (3 r is the only nonzero component 
of the shift vector and a is the lapse function. 

The procedure used to solve Einstein's equations is 
based on a 3+1 decomposition of the space-time and the 
solution of an initial value problem assuming a puncture 
type of initial data In order to follow the evolution 
of the initial data we found appropriate the use of the 
Generalized BSSN (GBSSN) evolution system of equa- 
tions defined previously for spherical symmetry in [5j as 
opposed to previous successful analyses related to the ac- 
cretion of a scalar field using Eddington-Finkelstein co- 
ordinates under the Arnowitt-Deser-Misner (ADM) for- 
mulation that requires excision Q. 

The evolution equations in the GBSSN system for the 
spherically symmetric space-time described by ([2]) reduce 
to a set of evolution equations for the conformal fac- 
tor Xi t ne conformal metric functions g rr and ggg, the 
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nonzero trace-free part of the conformal extrinsic curva- 
ture A rri the trace of the extrinsic curvature K and the 
contracted conformal Christoffcl nonzero symbol T r . The 
explicit expression for the evolution of these variables can 
be found in [1, 0] . 

The remaining equations of the geometry involve the 
evolution of the gauge. For the lapse we choose the 1+log 
slicing condition dta = /3 a d a a—2aK, and for the shift we 
implemented the recipe for the T-driver condition dtf3 a = 
lB a + (3 c d c [3 a and d t B a = d t T a +(3 c d c B a -f3 c d c r a -nB a , 
which helps to avoid the slice stretching near the hori- 
zon which is known to kill the numerical evolution. In 
order to avoid instabilities near the puncture we imple- 
ment a sort of excision without excision using a factor 
function on the sources of the evolution equations of the 
form (r/(l + r)) 4 from the coordinate origin out to one 
quarter of the size of the apparent horizon radius. Even 
though this function violates the constraints, the conver- 
gence tests show that such violation does not propagate 
outside of the black hole horizon. 

The evolution of the scalar field is given by the Klein- 
Gordon equation 



a<j) = 



-gg^d^} = -d+V, 



(3) 



where g is the determinant of the space-time metric. Wc 
solve this equation as a set of two equations for two first- 
order variables, tt = dt4> and £ = d r <p, coupled to the 
evolution of the geometry of the space-time. 

Initial data. In order to add the contribution of the 
scalar field to Einstein's equations, it is necessary to solve 
the constraints at the initial time slice which we assume 
for simplicity to be time-symmetric, which in turn im- 
plies that the momentum constraint is identically satis- 
fied. The scalar field contribution consists of a scalar field 
pulse which is used to solve the Hamiltonian constraint. 
Another important ingredient is that we set a space-time 
ansatz similar to that of a Schwarzschild black hole in 
isotropic coordinates. Thus we assume that the metric 
at initial time has the form 



fjr 



1, gee = r 
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(4) 



where M is the ADM mass at the puncture and agrees 
within numerical precision with the mass of the appar- 
ent horizon and u = u(r,t = 0) is the function to be 
determined through the solution of the Hamiltonian con- 
straint. The Hamiltonian constraint in terms of u reads 



d rr u = -(<9 r < 



2d r u 
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(5) 



is also possible to construct initial data consisting of a 
large enough amount of scalar field so that the ADM 
mass could become negative. 

It only remains to choose the gauge at initial time, 
for which we use an initially precollapscd lapse a = (1 + 
M/2r)~ 2 , zero shift (3 r = and zero time derivative B r = 
0. Because the initial time slice has been chosen to be 
time-symmetric, the extrinsic curvature components are 
all set to zero initially. These data are used to start up a 
Cauchy evolution using the GBSSN evolution equations, 
the gauge conditions and the Klein-Gordon equation. 

Numerical methods. The numerical method used to 
approximate the constraint and evolution equations is a 
second-order finite differences approximation. We only 
perform the evolution on a finite domain with artificial 
boundaries at a finite value of the radial coordinate. The 
integration in time uses a method of lines with a fourth- 
order accurate Rungc-Kutta time integrator. Through- 
out the evolution wc monitor the Hamiltonian, Momen- 
tum and Gamma constraints We check that they 
converge with second order to zero as the resolution is 
increased. 

Apparent horizon location. As we are interested in 
tracking the behavior of the black hole horizon, a useful 
diagnostics tool is the location of the apparent horizon 
during the evolution. We locate it from the definition of 
a marginally trapped surface (MTS) through the condi- 
tion = Viri 1 + Kijn % n' J — K = 0, where n l is an out- 
ward pointing unit vector normal to the apparent hori- 
zon, Kij are the components of the extrinsic curvature 
of the spacelike hypersurfacc on which one calculates the 
MTSs, and K its trace. The apparent horizon is the out- 
ermost MTS. For the metric ^ this equation is written 
as 
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(x(d r gee) - gee(d r x)) 
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(6) 



In order to track the evolution of the apparent horizon 
we calculate at every time step and locate the out- 
ermost zero of it at the coordinate radius rAH and cal- 
culate the mass of this horizon Mah — Rah/2, where 
Rah = v gee/x 1S the areal radius evaluated at tah- 

Mass. In order to have an estimate of the space-time 
mass we use the Misner- Sharp mass defined by @ 
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(7) 



which we solve using an ordinary differential equation 
integrator. We finally rescale r in such a way that x ~ * 1 
when r — > oo. We notice that with this procedure it 



where R = \J gee/x is th e areal radius. 

We estimate the ADM mass as Madm = 
linv^oo Mms so that we can estimate the mass of the 



3 



space-time at any time during the evolution even though 
we know it should be constant in time. 

Null rays and event horizon. We track a bundle of ra- 
dial outgoing null rays in order to i) approximately locate 
the event horizon as the surface that separates outgoing 
null rays that reach spatial infinity from those that do not 
when launched toward the future and ii) certify that the 
light-cones point in the correct direction at the apparent 
and event horizons. 

Results. We analyze different cases with various pa- 
rameters of the initial scalar field profile and different 
values of A including the case A — and obtain qualita- 
tively similar results. We use a Gaussian profile or a train 
of successive Gaussians with different amplitudes and 
widths for the scalar field. In particular we present two 
cases characterized by positive and negative ADM masses 
0. In both cases A = —0.1; in the case of Lagrangian (TT]) 
a negative potential with a maximum is appropriate to 
keep the scalar field stable. The results are shown in Fig. 
[1] where we show the evolution of the apparent horizon, 
Misner-Sharp and ADM masses 0. In both cases the 
area of the apparent horizon decreases and asymptoti- 
cally in time converges to the Misner-Sharp mass. It is 
worth noticing that for the case of negative ADM mass 
we used an initial train of scalar field Gaussians that are 
being successively accreted, which produces the step-like 
mass decrease of the horizon mass observed in Fig. 1. 

In order to show that the decrease of the black hole 
area is not a gauge artifact, in Fig. [5] we show a fine-tuned 
bundle of outgoing null rays for the two cases, and show 
the approximate location of the event horizon. The fact 
that the event horizon decreases with the accretion of the 
scalar field is a strong evidence that the hole area is truly 
decreasing since the event horizon is a gauge invariant 
surface unlike the solely consideration of the apparent 
horizon. 

Another interesting result is that the accretion rate 
and localization of the scalar field depends on the scalar 
field potential. As an example, in Fig. [3] we compare 
the behavior of the scalar field in the case where the 
potential is nonzero with the case of zero potential. In 
the case of the zero potential the scalar field is quickly 
accreted whereas in the other case the scalar field gets 
packed near the horizon and some scalar field is radiated 
away. 

Finally, in order to validate our results, we have verified 
the convergence of our results to second order so as the 
convergence of the violation of the constraints to zero. In 
Fig. [4] we show the convergence of the L2 norm of the 
constraint violations to zero for the simulations presented 
in Fig. [U 

Discussion. In order to illustrate the process we con- 
sider the hypothetical case of a black hole space-time with 
positive ADM mass. Assuming the mass of the original 
black hole apparent horizon is M p } lys — 10 9 M Q , and that 
a pulse of scalar field with amplitude 0.25 is accreted, the 
time-scale for the field to reduce the mass of the hole in 
25% is t ~ lday in coordinate time (see the results in 
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FIG. 1: We show the evolution of different masses versus 
the coordinate time. The solid line represents the mass of 
the apparent horizon, the dotted line the Misner-Sharp mass 
extracted at r ~ 15M% H and the dashed line shows the ADM 
mass of the space-time. In the left panel the ADM mass of 
the space-time is positive and in the right one it is negative. 
In both scenarios the apparent horizon mass decreases and 
approaches to the Misner-Sharp mass. In the first case the 
black hole loses approximately 25% of its mass while in the 
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FIG. 2: We present a bundle of outgoing null rays showing the 
approximate location of the event horizon for the runs shown 
in Fig. [T] It can be seen that this horizon decreases so as 
the apparent horizon does, showing that our results are not 
slice dependent. It is interesting to notice that the apparent 
horizon lies outside the event horizon. 

Fig. rj. 

In theory, due to the fact that our results scale with the 
initial mass of the black hole, it would always be possible 
to construct the appropriate scalar field profile so that 
the process of accretion could be iterated reducing half 
of the mass of the horizon each time. 

The results presented constitute solid steps toward ex- 
ploring the accretion of phantom scalar field dark energy 
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FIG. 3: We show snapshots of the scalar field at various times 
for A = (continuous line) and A = —0.1 (dashed line). The 
initial profile of the scalar field is the same, although the 
energy density is slightly different for both cases, nevertheless 
the time scales are comparable. We point out various effects: 
i) there is an initial outgoing pulse of scalar field due to the 
time-symmetry of the initial data in both cases, ii) however 
in the A = case the amplitude of the field quickly stabilizes 
around zero in the domain left behind by the outgoing pulse, 
whereas the field does not stabilize in the nonzero case, and 
instead new pulses of the scalar field are being ejected from 
the region near the horizon. 



under different astrophysical conditions, mainly those re- 
lated to cosmologically motivated potentials and wave- 
lengths of the scalar field ■ It would also be interesting 
to investigate the full nonlinear accretion of phantom en- 
ergy in black holes immersed in a FRW geometry like 
those in Mc Vittie-type black hole solutions [l(| , because 
these scenarios could have cosmologically interesting im- 
plications. 

One of the properties of this type of field in the cosmo- 
logical scenario is that if is the stress-energy tensor, 
the ratio between the radial pressure p r = n^T^ and the 
energy density p = v) 1 n l 'T iU/ measured by Eulerian ob- 
servers with the normal to the spatial hypersurfaces, 
provides an equation of state such that uj r = 2n < — 1 in 
the case of nonzero potential, and u) r = 1 in the case of 
zero potential. These properties change when the config- 
urations depend on the spatial coordinates and the po- 
tential is nonzero. For instance, in the case discussed 
here, the scalar field depends on the spatial coordinate 
and there are regions where the equation of state u> < — 1 
is satisfied, however there are also regions where <p a P~ 



proaches zero and the potential becomes zero too, in 
which case the scalar field would behave as stiff matter 

UJ=l. 

These results should imply important restrictions on 
the parameters of phantom dark energy models related 
to the abundance and mass of black holes and the abun- 
dance and lifetime of PBHs for instance. We also expect 
that the type of process shown in this article competes 
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FIG. 4: We show the convergence of the L2 norm of the vio- 
lation of the Hamiltonian, momentum and F r constraints to 
zero. The plots a,b and c correspond to the simulation with 
positive ADM mass and the plots d,e and f to the negative 
ADM mass case. 



with other horizon- reducing processes of black holes. 

A detailed study containing various initial condi- 
tions, more violent horizon-shrinking processes, dynami- 
cal properties of horizons and the study of the possibility 
of producing a naked the singularity is in progress [ill ]. 
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